Model studies of diffusion in glassy and polymer ion conductors 
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ABSTRACT 

Ion-conducting glasses and polymer systems show several characteristic pecu- 
^ liarities in their composition-dependent diffusion properties and in their dynamic 

^ . response. First we give a brief review of the current understanding of the ion dy- 

namics in network glasses in terms of stochastic theories. Secondly, a model for 
PEO(polyethylene-oxide)-based polymer electrolytes is described. The equation of 
state is calculated from the quasichemical approximation. Using this information 
together with Monte Carlo simulations of diffusion at constant volume, we obtain 
' constant -pressure results for polymer chain and ion diffusion as a function of ion 

concentration. This theory allows us to make comparison with experiments which 
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are carried out at constant pressure. 



1. INTRODUCTION 



a I Ion conducting glassy^"^ and polymeric^' ^ materials show a large variability in their chem- 

^ ■ ical composition, local structure, ion diffusion properties and in many other physical 

properties.^ Investigation of these amorphous systems and of the relevant ion transport 
^ ■ mechanisms thus constitutes a wide and active area of research. 

o ■ 

Generally, ion migration in these materials is connected with slow relaxation processes 
on time-scales much larger than microscopic times. A well-known example occurs in 
^ i glasses, whose low-frequency dynamic response appears to slow down continuously as 

5^ ■ temperature is lowered.^' ^ A quite different situation is known for polymer electrolytes, 

where ions are strongly coupled to the polymer network. Therefore, ion diffusion slows 
down and gets almost suppressed when approaching the glass transition temperature.^' 
Theories, which aim at a description of ionic motions on such extended time-scales, will 
have to eliminate all fast (e.g. vibrational) degrees of freedom from the outset. One then 
arrives at coarse-grained models based on stochastic ionic moves in an environment that 
reflects the characteristic local structure of the system. 

Adopting this point of view, we shall first focus on several transport anomalies in glasses 
and discuss the respective state of stochastic modelling (section 2). In section 3 we intro- 
duce a lattice model for polymer electrolytes which emphasises the mutual influence of 
ion and polymer chain dynamics. ^^'^^ Monte Carlo simulations of diffusion in an (NVT)- 
ensemble will be supplemented by calculations of the equation of state via the quasichem- 
ical approximation (QCA) in the version of Guggenheim^^ and Barker. Combination of 
both methods will enable us to obtain transport properties in an (NpT)-ensemble. Con- 
sideration of a constant pressure p is important when comparing theory with experimental 
trends under varying composition, such as the increase in the glass transition temperature 
or the decrease in ionic mobilities upon adding salt.^°'^^ 



2. ION TRANSPORT ANOMALIES IN GLASSES 



Alkali-doped network glasses (oxide- or sulphide-glasses, like the systems M2O — Si02 or 
M2S — SiS2, where M represents the alkah ion) exhibit several peculiar features in their 
diffusion properties: 

i) The activation energy E{c) for dc-transport markedly decreases with increasing al- 
kali concentration c for dilute systems (small c) , while it saturates for strongly doped 
samples. This behaviour results in a very sensitive c-dependence of the ambient- 
temperature dc-conductivity.^ 

ii) The Haven-ratio H{c) shows a steep decrease with c and also saturates for larger 

iii) Conductivity dispersion is an important effect where, apparently, different mech- 
anisms are active. At lower frequencies, Jonschers' 'universal ac-response' applies 
(power-laws in frequency, with exponents n„ ^ 0.6 to 0.7^'''^^), while at high fre- 
quencies one observes a 'nearly constant-loss'-type response, characterized by an 
exponent rii ~ 1.^'^^ When temperature gets lowered, the latter process dominates 
the conductivity spectra. 

iv) Dynamical processes which govern the ac-conductivity spectra also manifest them- 
selves in quasi-elastic neutron scattering, ultrasonic attenuation and nuclear-spin 
relaxation. 

v) Mixed-alkali systems, e.g. {Na20)x{K20)i-x^S'i02, show a highly nonlinear depen- 
dence of Na— and K-diffusion constants and of other physical properties on the 
ratio of mixing x.^^ 

Up to now, no theory exists that accounts for all of these features. Nevertheless, deeper 
insight into their microscopic origin has emerged from specific model studies. Here, wc 
only briefly address some of the most recent approaches and refer the interested reader 
to the original literature. 

Properties of dispersive transport and composition-dependent effects in single-alkali 
glasses can be understood with reasonable accuracy from stochastic lattice gas models 
based on structural disorder and Coulomb-interactions. ^^"^^ In particular, the 'counter- 
ion model' has proved successful in explaining the behaviour of weakly doped glasses. It 
also proposes a mechanism for the 'nearly constant loss' response in terms of correlated, 
dipolar reorientation processes. On the other hand, for explaining mixed-alkali-effects, 
a quite different approach in terms of the 'dynamic structure model' has been suggested, 
which involves the concept of lattice mismatch and a dynamical adjustment of the local 
structure to each ionic species.^ '''^^ The design of a unique model by a proper combination 
of these ideas remains a task of future research. 
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3. LATTICE MODEL OF POLYMER IONIC CONDUCTORS 



As pointed out in the Introduction, ionic and network degrees of freedom in polymer 
ionic conductors are strongly coupled and influence each other in a fundamental way. 
For example, in PEO (polyethylene-oxide)-based electrolytes, binding of cations to the 
electro-negative oxygen atoms in the chain molecules can introduce crosslinks into the 
polymer network. Local fluctuations of chains, which are essential for providing open 
pathways for the ions, therefore get reduced and ionic mobilities decrease. At the same 
time the network viscosity rises. Clearly, with increasing ion concentration c these effects 
get enhanced, but they are expected to saturate when the number of cations approaches 
the number of binding sites. Some of these aspects have been studied recently with the 
aid of stochastic models consisting of an interacting system of lattice chains and two 
species of point particles, representing cations and anions, respectively. Details of this 
model and the simulation technique are described in Ref..^^ Our aim here is to investigate 
a somewhat simplified model which contains only one species of point-particles ("ions") 
which can bind to the chains. This simplification, of course, ignores all features connected 
with anion diffusion, yet it allows us to investigate the profound influence of ions on the 
polymer network rigidity. In addition, it greatly facilitates the problem to convert results 
for diffusion at constant volume to the corresponding results for constant pressure. This 
will enable us to describe experiments on polymer electrolytes under varying salt content, 
which normally arc carried out under constant pressure. 

The model we consider involves lattice polymer chains (on a simple cubic lattice of spacing 
a) which consist of two types of beads, C-beads which do not interact with the point- 
particles and X-beads which attract point-particles with strength e. In analogy to PEO 
the sequence of beads is taken as C(XCC)„, so that the length of chains is given by 
r = 3n + 1. Throughout this work we assume r = 13. All the beads repel each other 
with a common strength, which again is taken as e. In our simulations, sec section 5, the 
equilibration procedure and the rules for elementary chain and ion motions are taken as 
in Ref..^^ 

4. EQUATION OF STATE 

To extract the pressure from Monte Carlo runs for a lattice system normally is a difficult 
task, although several methods are known in principle. One reason is the problem 
of thermahzation following a stochastic volume change. Therefore, with respect to the 
equation of state of our model, we shall rely on the quasi-chemical approximation (QCA), 
which requires solving a set of nonlinear equations for the quantities N^j , which denote the 
number of nearest neighbour pairs consisting of molecules of type i and j. The version by 
Barker, which we adopt here, allows us to deal with heterogeneous chains as described 
in the foregoing section. 

Rather than going into details, we merely give the final result for the pressure p, obtained 
from 
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Here M, Nq and A^oo denote the numbers of lattice sites, vacant sites and vacancy pairs, 
respectively, to be calculated from the QCA-equations as a function of the numbers of 
molecules and temperature. 

Isobars for different numbers of ions are shown in Fig. 1 at a fixed (reduced) pressure 
pa? /e = 0.35. When we add ions at a high temperature, the system swells, but the amount 
of swelling is less than what one would expect for non-interacting point particles. This, 
clearly, is a consequence of the attractive ion-polymer interaction, which even leads to a 
reversed trend at low temperatures (sec Fig. 1). Let us remark that inclusion of anions in 
an extended model would favour a volume increase. 
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Figure 1: QCA-isobars {pa^/e = 0.35) for different concentrations of point particles. Values for volume 
and temperature indicated by special points are used in the simulations for and D'^^^ shown in Fig. 2. 
At high temperatures (e/fcsT < 1) the different curves refer to (from below): A^i = 0, 20, 40, 60, 120, 240. 

5. DIFFUSION PROPERTIES UNDER CONSTANT PRESSURE 

Next we turn to Monte Carlo simulations of the long-time diffusion constants D^^^ and 
L>(^) for ions and for the centres-of-mass of polymer chains, respectively. These quantities 
are deduced as usual from the corresponding time-dependent mean-square displacements. 



Thereby we choose a system with Np — 31 chains of length r = 13 and with a varying 
number A'^i of ions. Given a fixed pressure p and temperature T, we carry out these 
simulations by choosing a nearly cubic simulation box of size Li x L2 x L3 = V{p,T)/a^ 
(with periodic boundary conditions) such that the volume V{p, T) corresponds to that 
obtained from the QCA (see Fig. 1). Fig. 2 shows the T-dependence of both D'^^^ and D^^^ 
at a pressure given by pa^ /e = 0.35 for various amounts of ions. As T goes to infinity, we 
approach the dilute limit {V 00), and chain diffusion constants for different Ni assume 
a common value (that corresponds to the athermal limit). As temperature is lowered, the 
different curves in Fig. 2 show a downward curvature of a form which is well represented 
by the Vogel-Tammann-Fulcher (VTF) law, 

10° I ^ , ^ , 1 




Figure 2: Simulated polymer chain and ion diffusion constants, D*^^) and D^-^\ versus inverse temperature 
at a constant pressure, pa^/e = 0.35. The different sets of data and (fitted) curves refer to different 
numbers of ions (from above): (iVi = 0(0), 20(A), 40(v), 60(t>), 120(0), 240(v). At a given temperature, 
the sizes of simulation cells were chosen according to the volumes deduced from the QCA, see the special 
points and corresponding symbols in Fig. 1. Continuous lines represent fits in terms of the VTF-law. 



D'«'(r,iVO = fl£)exp(-j,^llilL_) (2) 

Here a stands for / and P. The quantities D^'' denote the diffusion constants at infinite 
temperature, D^^ = Dq being the diffusion constant of a free particle, and /Dq ~ 
10^^. £'"(iVi) is an energetic parameter and T^^'^Ni) the VTF-temperature. 
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A first important observation is that both D^^^ and D^^^ are reduced when wc increase 
A^i, as a result of an increased number of cations which bind to the chains (X-beads) and 
possibly form crosslinks. As a consequence, we observe rising VTF-temperatures T*^")(iVi), 
as displayed in Fig. 3. Generally, the VTF-temperature represents a lower bound to the 
glass transition temperature Tg. Our simulations with respect to the diffusion of chains 
therefore are consistent with the experimentally observed increase in Tg with the ion 
content. Interestingly, for weakly doped polymers (A^i^60), T^^\Ni) and T^^\Ni) roughly 
coincide, which is indicative of a strong coupling between ions and chains. However, as A^i 
becomes comparable with the number of X-beads, which in the present case is ANp — 124, 
the influence of Ni on the ion and chain mobilities diminishes, as seen by comparing the 
data for Ni = 120 and Ni = 240 in Fig. 2. Simultaneously, the VTF-temperatures 
saturate, as shown in Fig. 3. At saturation, T^^^ < T'^^\ which means that for T ~ T^^^ 
the ions maintain a certain mobility while the system of chains freezes. The reason is that 
in the strongly doped regime only part of the ions will bind to the X-beads. These results 
correctly describe the experimental trend in the glass transition temperature of PEO- 
based electrolytes as a function of ion content in an appropriate regime of doping. -"^^'^^ 
At even larger ion concentrations, however, the experimental Tg drops after passing a 
maximum. This is possibly due to ion clustering effects, not taken into account in the 
present model. 
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Figure 3: VTF-temperatures for ions (T(^)) and polymer chains (T^^)) versus the number of ions, as 
deduced from the data in Fig. 2. (Continuous Hues are guides to the eye.) 
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6. SUMMARY AND OUTLOOK 



Diffusion properties of amorplious solid electrolytes were discussed, emphasizing the ac- 
tual capabihties of semi-microscopic (stochastic) models and Monte Carlo simulation to 
account for the experiments. Naturally, in this context we were interested only in experi- 
mental key features, which are common to many different materials. 

After pointing out several characteristic ion transport properties of glasses together with 
successes and persisting problems in their theoretical understanding, we turned to a spe- 
cific model for PEO-based polymer electrolytes. Composition-dependent diffusion proper- 
ties and the tendency to undergo a glass-transition were studied under constant pressure. 
This was achieved by combining constant-volume simulations with the equation of state 
derived from the QCA. 

The model we studied contains only one species of particles (cations) besides the polymer 
chains and hence does not allow us to investigate simultaneously the diffusion of anions. 
A simplified treatment of anion diffusion in an otherwise realistic model of a polymer elec- 
trolyte might be developed with the help of the concept of dynamic bond percolation. ^^"^^ 
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